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SUMMARY 


Flight tests were conducted for the purpose of verifying an analytical 
aerodynamic derivative model of a CH-47 tandem-rotor helicopter at low cruise 
speeds and transition to hover portions of curved, decelerating landing approach 
flight paths. The flight testing was performed on a closed loop basis with the 
stability augmentation system (SAS) of the CH-47 operating, and transient 
response data were obtained using both manual and computer -generated input 
maneuvers. 

The model verification consisted mainly of comparing the existing analyti- 
cal derivatives with those identified from the flight data using both extended 
Kalman filter (EKF) and maximum likelihood estimation (MLE) algorithms, then 
comparing the measured response time histories with those predicted by each of 
the three sets of derivatives. The overall evaluation further included verifi- 
cation of the analytical SAS model, as well as investigation of the effects of 
closed loop flight testing on derivative identification. 

The results indicate some amplitude and frequency differences between the 
measured response time histories and those predicted by the analytical deriva- 
tives that vary in magnitude with each test run; these differences appear to 
be mainly due to inaccurate values for some derivatives. The presence of non- 
zero trim accelerations in the flight data affected the MLE and EKF identified 
derivatives so that the time histories predicted from them also differed from 
the measured ones. With some exceptions the discrepancies are not severe, and 
the overall agreement between the measured and computed time histories is rea- 
sonably good. 

The results further indicate no adverse effects attributable to closed loop 
flight testing, nor any deficiencies in the analytical modeling of the SAS. 

The use of computer-generated input maneuvers proved to be superior to manual 
ones and was found to be highly effective in generating transient response data 
having good modal excitation without exceeding small perturbation amplitude 
bounds. Subsequent flight testing for other purposes indicated that the prob- 
lem with nonzero trim accelerations could be largely eliminated by also using 
the fly-by-wire control system to automatically trim the aircraft prior to 
initiating a test maneuver. 


INTRODUCTION 

Avionics research for helicopters has been in progress at the Langley 
Research Center over the past several years as part of the VTOL approach and 
landing technology (VALT) program (ref. 1). An NAS A/Army /Boeing Vertol CH-47 
helicopter served as the basic research aircraft for a comprehensive test pro- 
gram which has as its ultimate goal the development of technology necessary for 
optimum VTOL short haul transportation in the coming decade. The need for a 
suitable analytical description of the helicopter dynamics, and the rationale 



for selecting a linear perturbation model based on aerodynamic stability and 
control derivatives, is delineated in reference 2. 

The extensive tables of aerodynamic derivatives contained in reference 2 
were generated by means of an analytical small perturbation technique, involving 
detailed equations for the helicopter forces and moments, and cover essentially 
the entire CH-47 flight regime. Although these derivatives appeared to provide 
a good representation of the helicopter dynamics, comparison with results 
obtained from actual flight data remained as an important step in verifying 
the chosen analytical model. The flight tests described in reference 3 were 
designed to obtain suitable data for this purpose and, also, to permit identi- 
fication of aerodynamic derivatives for comparison with the tabulated analytical 
values. Because landing approach and transition to hover from low cruise speeds 
of about 60 knots comprise the flight regime of primary interest for the VALT 
program, flight data for the remainder of the CH-47 operating range were unnec- 
essary and were not obtained. 

The actual flight testing and ensuing derivative identification activities 
were planned and conducted in cooperation with the Structures Laboratory of the 
U.S. Army Research and Technology Laboratories (USARTL) at the Langley Research 
Center. A joint effort was undertaken with the objectives being to make the model 
verification as comprehensive as possible and to provide a basis for evaluating 
differing methods of identifying aerodynamic derivatives from flight data. Ref- 
erence 4 contains comparisons of analytical derivatives from reference 2 with 
those computed by USARTL using an extended Kalman filter (EKF) algorithm and 
with preliminary values from the maximum likelihood estimator (MLE) employed 
in the present report. 

The primary purpose herein is to document the completed MLE results and 
to include comparisons of measured transient response time histories with those 
predicted by each available set of analytical and identified derivative values. 

The effect of the CH-47 stability augmentation system (SAS) on the derivative 
identification process was also investigated, as the flight testing was neces- 
sarily conducted on a closed loop basis to avoid unstable responses which 
quickly exceed acceptable small perturbation bounds. This mode of testing fur- 
ther provided the opportunity to verify the analytical model of the SAS given in 
reference 5. An additional objective of the flight test program was to compare 
results obtained by using conventional doublet type manual input maneuvers with 
those from computer-generated inputs, implemented by utilizing the fly-by-wire 
control system of the test aircraft in conjunction with an onboard computer. 

(See ref. 4.) 


SYMBOLS 

A,B,C,K,Z aircraft SAS model matrices (eqs. (10) and (11)) 
A-| ,A 2 control input amplitude coefficients, cm (eq. (21)) 
b measurement bias vector 

d data sample (eq. (23)) 
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F,G 


aircraft dynamics system matrices (eq. (4)) 


f frequency, Hz 

g acceleration due to gravity, m/sec^ 

H,D aircraft response measurement system matrices (eq. (3)) 

h n digital filter weights (eq. (26)) 

I X ,I Z ,I XZ aircraft moments of inertia, kg m^ 

J MLE performance function (eq. (1)) 

N number of data samples or digital filter weights 

n x , n y,n z body axis components of aircraft linear acceleration, g units 

P,Q,R roll, pitch, and yaw rates about aircraft body axes, respectively, 

deg/sec 

p ,q,r perturbations in P, Q, and R, respectively, deg/sec 

p MLE parameter identification vector (eq. (2)) 

R (t) weighting matrix 

S number of filtered data samples 

t time, sec 

U,V,W body axis components of aircraft linear velocity, m/sec 

u, v, w perturbations in U, V, and W, respectively, m/sec 

U aircraft control input vector 

u perturbation in U (eq. (10)) 

v aircraft SAS model state vector (eq. (10)) 

w white measurement noise vector (see equation following eq. (5)) 

X aircraft state vector 

x perturbation in X (eq. (10)) 

Y aircraft output response vector 

y perturbation in Y (fig. 2) 

Slonq longitudinal control input, cm 
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Scoll collective control input, cm 

^lat lateral cyclic control input, cm 

^rud differential control input, cm 

s 4> , SAS ' 5 e , SAS ' 5 4 j, SAS roll, pitch, and yaw SAS feedback command signals, 

respectively, cm 

k control input amplitude scale factor (eq. (21)) 

$,© roll and pitch attitudes about aircraft body axes, deg 

4>,9 perturbations in $ and 0, deg 

tp aircraft stability and control derivatives vector 

w frequency, rad/sec 

Subscripts: 

c control input or cut-off frequency 

f final value 

m measured value 

i,j,k indices for scalar, vector, or matrix quantities 
n index 

o initial, trim, or central value 

t termination frequency 

SAS stability augmentation system 

Notation: 

E ( ) statistical expectation 

( ) matrix inversion 
( ) T matrix transposition 

(^) estimated value 

A( ) increment 

(“) normalized or average value 
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Abbreviations: 

EKF extended Kalman filter 

MLE maximum likelihood estimator 

SAS stability augmentation system 

US ART L U.S. Army Research and Technology Laboratories 

VALT VTOL approach and landing technology 

VTOL vertical take-off and landing 


MAXIMUM LIKELIHOOD ESTIMATOR 

The aerodynamic derivative estimates documented in this report were com- 
puted by means of the maximum likelihood estimator of reference 6. The mathe- 
matical development of the MLE algorithm as used herein is described by equa- 
tions (1) to (9); and the aerodynamic derivative model employed for the CH-47, 
by equations (10) to (20). 


Algorithm Formulation 

The algorithm for implementing the MLE estimator is formulated in the same 
manner as in reference 7, except for added computational options, and is based 
on minimizing the error in the perturbation response output of the aircraft in 
the least squares sense by using a performance function of the form 


J(p) = ' y(P-.t)] T R(t) 1 [y m (t) - y(p,t)]dt 

to 


where y m (t) and y(p,t) are the measured and estimated output response vec- 
tors over a data-gather ing period [t 0 ,tf]], R(t) is a symmetric, positive- 

definite weighting matrix. Since the algorithm is developed in discrete form 
for computational reasons, J(p) is approximated: 


N 


j(p) - ^ [y m ,i - y(p>i] R 1 [y m , i - y(p)J 


(i) 


i=l 
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The corre- 


which represents N samples of the output error during 

3J 



sponding least squares normal equations — =0 are then solved for p, which 

8p 

includes the derivative estimates $, by means of the differential correction 
procedure 


Pj+1 = Pj + 4>j 


= Pj + 



( 2 ) 


which is also called a quasilinearization or modified Newton-Raphson minimiza- 
tion technique. (See ref. 6.) As in reference 7, the convergence criteria used 
was j Ap j j S | 0. 01 pj j simultaneously for all elements of p. 

If the measured control input and output response vectors u m ^ and y m ^ 
are assumed to have the forms 


-“m, i 


u i + ^c + w c, i 


and 


^m, i 


= y t + b + wj_ 


which account for the presence of biases and random white noise, the estimated 
output y^ required in evaluating equation (2) is then modeled as 


yi = H (p) + D (p) u t + b (3) 


where the state vector Xi is obtained by numerical integration of the aircraft 
dynamical equations 



F(p) Xi 


+ G(p) Ui 



(4) 
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A 

and the control input vector is estimated from 


A 


u 


i 



( 5 ) 


The weighting matrix R in equation (1) is conventionally taken to have the 
diagonal form 


R6j_j E[w i Wj] T 


where eQw^] =0 at every sample point. Furthermore, the estimate R used in 
equation (2) is recalculated at each iteration of the differential correction 
process from the expression 


'1 



N 



i=1 



(j = k) 

(j t k) 


( 6 ) 


Other than the specific forms of the vectors and matrices, which are defined 
with the aerodynamic derivative model for the CH-47 in the next section, the 
remaining quantities to be determined in evaluating equation (2) are the 
partial derivatives 3yf/3p obtained by differentiating equations (3) 
and (5) : 


9 Yi 3x i 0H(p) - 3D(p) 

= H(p) + X; + 

3p 3p 3p 3p 


D(p) 


3b, 

3b, 


c 3b 

+ — 

3b 


(7) 


in which 3b c /3b c and 3b/3b (as well as 3x/3x 0 in eq. (8) below) are equiv- 
alent to identity matrices of suitably corresponding dimensions. Since the full 


rank of 

ing the 
biases , 


the parameter vector p = ! fe c|bjx 0 jT includes the option for estimat- 

initial state conditions x 0 as well as the aerodynamic derivatives and 
the following partitioning of 


r\ ^ *1 A | /\ | A 

9 Y i 9y j 9y 1 3y J 3y 

3p 3cp I 3b c | 3b | 3x 0 
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and 


3 x i 9x| 9x J 9x 

9p 3cp|9b c |9x 0 


serves to indicate the corresponding component structure of these matrices. 
The values for 9x^/9p are generated by numerical integration of 



9x j 


3G(p) 


9 b, 


' x o * ( x ) t=0 


F(p) + x i + , Ui - G <P). k 

dp dp dp db c 

i 

9 x o 

t 1 A 

i i 9 X 


V 3 ? 

°! °| — 
u 1 1 3x °J 


( 8 ) 


t=0 


which is obtained by differentiation of equations (4) and (5) . 
9Yi/3p and 9x^/9p form according to the rule for Jacobians, 

9H(p) a 

for the matrices resulting from the tensor products x^, 

9F(p) , 9G(p) a 3p 

xj , and u.; are defined by 

9p 1 9p 1 


The elements of 
whereas those 


9D(p) a 



a ik 


9C 

— X 

9p 


I 


j=1 


9 c i j 

x i 

3 Pk 3 


(9) 


The algorithm is further configured to provide options for deleting measurements 
from y m ^ and for holding any combination of elements of p constant during 
the iterative solution process. To delete, say, the kth measurement from y m 
the kth row of 3y^/9p and kth element of (y ra ^ - y-jj are set to zero. Holding 
the rth element of P constant is accomplished in a similar manner by setting 
the rth row of 9y^/9p to zero, but with the additional requirement that the 
rth element of the rth row of the resulting information matrix 

A ' 

9 Yi 

9p 


5-1 


3y i\ 
>3P 


then be reset to have the value of unity instead of zero. 
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Aerodynamic Derivative Model 


The linear perturbation model of the CH-47 dynamics employed in reference 2 
is characterized by equation (4) , which expresses the aircraft response relative 
to an equilibrium or nominal trim flight condition in terms of the control input 
and state vectors. Equations (10) and (11) define the extension of equa- 
tions (3) and (4) to include modeling of the helicopter stability augmentation 
system (SAS) and are given in closed loop form by 



(Dju + b 


( 10 ) 


and 


d 

X 


; 

GC 

dt 

__ V 


_ KB | 

A 


+ [G]u 


f x o = (X) t=0 

\V 0 = (v ) t=0 


( 11 ) 


in which the rigid body and SAS portions of the model are indicated by partition- 
ing. The corresponding segments of the total state vector are denoted by x 
and v. The aerodynamic stability and control derivatives constitute the prin- 
cipal elements of the matrices H, D, F, and G, and those elements to be 
identified by means of equation (2) comprise the portion of the parameter 

vector p. The remaining matrix components GZ, GC, KB, and A pertain to 
the SAS model. 

While originally included to permit verification of the model obtained by 
converting the SAS transfer functions given in reference 5 to state variable 
form (see ref. 3) , the resulting SAS model also provides the useful option of 
performing derivative identification with either open or closed loop formula- 
tions. Equations (10) and (11) revert to open loop form by simply deleting the 
SAS terms; however, the effect of the SAS still must be taken into account since 
all flight tests were conducted on a closed loop basis with the SAS operating. 
This requirement is easily satisfied by adding the SAS output signals, which 
were measured as part of the CH-47 flight test data, to the control stick inputs. 
For verification against their measured values, analytical estimates of these 
commands are obtained from the SAS model by means of 


^9, SAS ~ V 1 

(12) 

5 <1>,sas = -K pP 

(1 3) 
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5 ^,SAS = K v v - K r r + v 5 + v 7 


( 14 ) 


where the functions of the gains and state variables appearing in these formulas 
are subsequently indicated in the appropriate ones of equations (15) to (18). 

The numerical values used for all gains and other SAS model constants are listed 
in table I. 

Although the derivatives tabulated in reference 2 are for the usual six 
degree-of-freedom formulation, the following expressions for the corresponding 
uncoupled longitudinal and lateral three degr ee-of-f reedom representations were 
used in generating the MLE derivative values presented in this report. This 
simplification was made on the basis of the conclusion that both simulation 
studies (see ref. 3) and the actual flight data indicated the SAS uncouples 
these modes very well. 

Longitudinal mode .- For the longitudinal motions of the helicopter, the 
respective three degree-of-freedom forms of equations (10) and (11) are 


0 


1 

0 

0 

0 

1 » 

0 

0 

q 


0 

1 

0 

0 

i • 

0 

0 

gn x 

= 

0 

x q 

x w 

x u 

1 

1 X ®long 

0 

0 

Wz 


0 

Z q 

z w 

z u 

1 2 ^long 

0 

0 



0 0 



b 0 


0 0 


^long 


b q 

+ 

X( ^long X< $coll 


^coll 

+ 

b n,x 


Z| ^long 2 ^coll 



b n, z 


v 3 


(15) 


and 


e" 


0 

1 

0 

0 

1 

1 

0 

0 

0 “ 


0 

q 


0 

Mg 

«w 

»u 

1 

1 

1 

M ^long 

0 

0 


q 

W 


-g sin 0 O 

Zq + °o 

Zw 

Zu 

1 

1 

Z ^long 

0 

0 


W 

U 


-g cos 0 O 

Xq ” W 0 

X w 

X U 

1 

1 

X ^long 

0 

0 


u 

VI 


0 

Kqbl 

0 

0 

1 

1 

0 

1 

0 


V 1 

v 2 


0 

Xqt>2 

0 

0 

1 

1 

0 

0 

1 


v 2 

_ v 3 __ 


0 

Kq b 3 

0 

0 

I 

1 

_a o 

- a l 

-a 2 


v 3 _ 


0 

0 

M ^long 

M|5 coll 

2 ® long 

Z| ^coll 

X ®long 

X ®coll 


^long 
_^co 1 1_ 


( 16 ) 


for which the individual elements of the parameter vector p components are 
defined by 
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9 - [Hq l WV|ZqZwZuXqV l u M <5 long M <5 co ll Z 5 long Z 5 coll X(5 long X(S coll] T 
b = [b0bgb nf x b n , z] T 

x Q = [eqwu ] t=0 

Lateral-directional mode .- The corresponding three degree-of-f reedom equa- 
tions for the lateral helicopter motions are 



where the individual elements of the p vector components are given by 



b <= - [ b «lat b 5 rucl] T 
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b = [bpbj-b^yJT 
T 

x o = [vpr<|>]t =0 

and 


Y v 

= Y v + K v Y 6rud ; 


Y P K P Yfi lat ; 

Y r = Y r - K r Y 6rud 


= L v + K v L 5rud ; 


= L p - K p L 6i at ; 

L*’ = L* - K r L^ ud 

< 

= < + k v n* w 

■S' 

= - K p N 6 lat ; 

N*' = N* - K r N5 rud 


which contain feedback terms corresponding to the respective gain elements in 
equations (13) and (14) for the roll SAS and yaw SAS commands. The roll and 
yaw derivatives are further modified according to 





1 

I xz 




1 


H 

L i 



I 2 

• L xz 

_ I X I Z_ 


J XZ 

J-z 

i 

_ N i 


(i = v,p,r, 5 lat , 6 rud ) (19) 


which account for cross products of inertia ratios that couple the p and r 
equations. (See ref. 8.) Conversion of the identified derivatives to the for- 
mat of reference 2 is given by the inverse relationships 


L 


i 





(i v , p, r , Slat > ^r ud ) (20) 


and the values given in reference 2 for the CH-47 moments of inertia are 
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Ix 


= 50 386.3 kg m 2 


I z = 257 685 kg m 2 


I X Z = 19 838.3 kg m 2 


Lastly, the trim values of U Q , W Q , and 0 O appearing in equations (16) 
and (18) are listed in tables II to IV for each of the flight test runs. 


PROCEDURE 

The CH-47 flight tests were conducted for the conditions given in table V, 
which are representative of the landing approach trajectories of interest for 
the VALT program. The 15 principal test points listed involve a total of 
44 individual runs and encompass most of the test plan presented in reference 3. 
As conventional pitot-static airspeed indicators are too inaccurate below about 
50 knots, the time histories of the helicopter velocity for each run were 
obtained instead by means of precision radar tracking. (See ref. 4.) The cor- 
responding wind magnitudes for correcting the velocity data to local freestream 
values were determined in the same manner by tracking weather balloons between 
successive runs and are also listed in table V. 


Input Maneuvers 

The test program set forth in table V includes runs for both manual and 
computer-generated input maneuvers. Due to the complexity of the piloting task, 
the former were limited to single control stick commands while the latter were 
implemented with both single and double ones as indicated. The manual inputs 
are of the usual doublet type, whereas those generated by means of the onboard 
computer consist of a long and a short period sinusoid described by 


u(t) = K (A-] sin o>i t + A 2 sin c^t) 


( 21 ) 


and are depicted by the sketch accompanying table VI which lists the amplitude 
and frequency coefficients for each of the four basic helicopter control stick 
commands 5 long' 5 co ;li, ^lat' and ^rud* T ^ e frequencies of these sinusoids, 
for both the longitudinal and lateral dynamical modes, were chosen to correspond 
to the principal long and short period characteristic roots of equations (16) 
and (18). Each of the single input maneuvers for the runs at hover and at 
20 knots were generated by evaluating equation (21) for only one of the four 
basic control stick commands, while the double input ones for the 40 and 
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60 knot runs employ a second such sinusoidal stick command commencing 6 sec 
later than the first one. The values entered in table VI for the starting 
times t 0/ i and t 0r 2 permit combining ^coll with S^ong anc3 ^rud with 
Slat as required for the runs listed in part (b) of table V. To avoid nonlin- 
ear responses, A 2 sin u^t of ^long for the double input maneuvers is delayed 
6 sec so that the two sinusoids occur separately instead of simultaneously. 

In addition to the scale factor K, which was included with the computer 
implementation of equation (21) to permit varying the amplitudes of the u(t) 
so as to restrict the y m , i to linear values, a pilot stabilization model of 
the form 


5 longpiiot = K<°- 2117 u " 0.7867w) 


( 22 ) 


was also programed. The use of 


6 longpiiot 


was recommended to counter pitch 


response to S-L at inputs, predicted by simulation results described in refer- 
ence 3, that can cause large nonlinear changes in u of 9 m/sec or more at 
airspeeds of 40 knots and higher. According to reference 3, this response will 
occur because the long period pitch divergence mode becomes appreciably unstable 
at airspeeds of about 40 knots. 


Data Processing 

After completion of the basic data reduction, the flight data records were 
edited and then filtered to remove high frequency noise prior to arranging them 
in perturbation form for use in performing the derivative identification and 
model verification tasks. The filtering was accomplished by means of the digi- 
tal filter of reference 9 


N 


df - h Q d^ + ^ h n (di._ n + ^i+n) 
i=l 


i > n 
i = 1 to S, 


(23) 


in which 


^i-n " ^i+n 


1 s n 

i = 1 to S 


i > (S - n) 


^i+n ^i-n 


1 4 


i = 1 to S 


and and d^ denote the filtered and unfiltered values of the ith sample 

of a measured quantity in a data record containing a total of S samples. 
Equation (23) is an alternate formula for the filter that is convenient for 
handling the beginning and ending portions of the data record where values of 
d^ that do not lie within Q,s3 are called for in the filtering process. 

The weights h k apply to 2N + 1 of the S samples for each value of i 
and are normalized according to 


h k 

h k = (24) 

N 

ho + 2 J h n 

n=l 


in which each weight is divided by the sum of all 2N + 1 values. The number 
of weights to be used, which depends on the cut-off and termination frequencies 
f c and f t and the sampling interval At, was determined by 


2 

N = 

(f t “ f c ) At 


(25) 


as recommended in reference 9. 

The set of weighting functions h n are generated by means of the formula 


1 


h n = 


2 Tin At 


sin ( 0 t n At + sin « c n At 
1 - 4(f t - f c ) 2 (n At) 2 


(n = 1 to N) 


(26) 


where w = 2iTf. As equation (26) is indeterminate for n = 0, and for 
n At = — - as well, the expression 


2 ( f t - f c > 


h n = 


ft- cos 'jJ t n it + f c cos 0 J C n At 
1 - 1 2(f t - f c ) 2 (n At) 2 


(n = 0 to N) 


(27) 


obtained by means of L'Hopital’s rule is then used. Evaluation of equa- 
tion (27) for n = 0 yields h Q = f fc + f c for the value of the central weight. 



A typical example of the filter effectiveness in removing high frequency 
noise is illustrated by the solid curve plotted in figure 1 , which represents 
the filtered roll rate time history from run 10. (See table IV.) The sampling 
rate and the cut-off and termination frequencies used to produce d^ for this 
plot, and for all of the flight data, were At = 1/40 sec, f c = 1 Hz, and 
ft = 2 Hz. These choices resulted in the use of 161 weights since the corre- 
sponding value of N, as determined from equation (25) , is 80. 

After applying the wind corrections from table V and resolving the result- 
ing freestream velocity components along the aircraft body axes, the remaining 
processing step was to represent the data as perturbations relative to equilib- 
rium or trim flight conditions. This reduction only amounts to forming the 
differences between the total output response measurements Y m ^ and their 
respective trim values Y mf 0 as indicated in figure 2. 

For the ideal situation where the choice of Y m 0 is free from error, 
the desired output response perturbation data is simply y mf i _ Y m, i Y m, o 
which also subtracts out any measurement bias. However, if problems such as 
unsteady trim conditions cause the selection of the trim point value to be in 
error as denoted by Yj^ 0 in figure 2, then Y m,o " Y m,o b and the resulting 
perturbation data become 


Y m, i = Y m, i 
= Y m, i 
= Y m, i 

so that the effect of the trim point selection error is to introduce a bias b 
into all i values of y^ i- Note that the small At 0 time offset between 
YjJj 0 and Y m 0 has no effect other than to shift the time at which y^, ^ 
starts and does not contribute any lead or lag. The bias estimate term in equa- 
tion (3) was included mainly as a means to compensate for the trim error biases 
just described. Equation (5) was similarly structured since essentially the 
same situation occurs with the control input measurements. 


■■m, o 


Y m, o + b 


+ b 


MI.E Derivative Identification 

The computation of 9 by means of the MLE algorithm given by equation (2) 
was not entirely successful since sane of the test runs in table V were compro- 
mised by dropouts in the radar measurements and failure to achieve good trim 
conditions prior to initiating the test maneuvers. 

The dropouts often caused uncertainty in determining the trim values of 
U c and W 0 required in evaluating equations (16) and (18) and precluded 
obtaining any values for them at all for some of the runs. This problem further 
prevented including the velocity perturbation measurements u^ = - U D , 

v^ = Vj^ - V Q , and w. = W-[ - W Q as elements of the y mf i vector in the appro- 
priate one of equations (15) and (17). Inspection of equations (15) to (18), 
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however, shows the derivative information contained in these measurements essen- 
tially duplicates that of the linear accelerations n x , ny, and n z so that 
omission of u, v, and w does not impair or prevent derivative identifica- 
tion. Thus, no further attempt was made to use u, v, or w as response mea- 
surements in generating any of the MLE results. 

The consequences of inadequate trim conditions gave rise to a more serious 
problem which in effect is equivalent to modeling error. Unsteady trims can 
introduce appreciable linear and angular acceleration contributions to the air- 
craft response motions that are assumed to be negligible in postulating the 
linear perturbation model. (See ref. 8.) Since the dynamical models defined 
by equations (16) and (18) do not correctly describe the aircraft response under 
such conditions, the MLE algorithm then cannot yield accurate derivative esti- 
mates. The trim conditions assumed in stating these equations are that 
V Q = P Q = Q 0 = Rq = 0 and that U 0 and W Q have constant values. Typical 
values for these quantities are illustrated by the curves in figure 3 which are 
plotted for run 18 from test point 5. (See table V.) 

Due to the changing values of the trim variables, the magnitudes of the 
corresponding acceleration contributions can vary sharply and can cause the 
convergence characteristics of ip to be very sensitive to the choice of the 
initial time t 0 from which the test maneuver is presumed to start. A change 
in t Q of only 1/40 sec, for example, caused divergent solutions to converge 
and vice versa for some of the test runs. The MLE derivative identification 
computations with the present flight data usually yielded convergent solutions 
for tp when t Q could be chosen such that the trim values of the linear and 
angular accelerations were within about 0.3 m/sec 2 and 0.5 deg/sec 2 , respec- 
tively. The slopes of the curves plotted in figure 3 show that these limits 
cannot be satisfied for any usable choice of t Q , which apparently is the pri- 
mary reason a convergent solution was not obtained for run 18. The best (p 

attempt occurred for t Q = 2 sec as marked on the plots in figure 3 but 

diverged on the 20th iteration. 

Unsteady trim conditions also appeared to be causally related to the slow 
convergence observed for several of the other test runs, which is characteristic 
of the solution behavior in the presence of modeling error. Convergence of 
equation (2) typically required 10 to 20 iterations, compared with 5 to 8 itera- 
tions for solutions based on the use of simulated y m i data which contained no 
modeling error. Coincidentally, the values of <p were often distorted from rea- 
listic levels, and those for £, b c , and x 0 were usually much larger than the 
maximum error that could possibly be introduced in forming y m i and u m 

The fact that including £ and b c in p caused the solutions for some runs 

to converge - albeit to physically incorrect values which were nonconvergent 
when only the cp were estimated - further indicates modeling error. Addition- 
ally, the values for <p, b, b c , and x Q were all correctly and accurately 
determined when identified from simulated data containing only measurement 
error . 

These results are consistent with the well known fact that the minimiza- 
tion process embodied in equation (2) will vary all of the elements of p, in 
accordance with the analytical structure of the x^ and y^ models, so as to 
minimize (y m , i - y^) in the least squares sense. (See ref. 10.) When modeling 
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error is present, the estimated parameters therefore become more and more dis- 
torted as the rank of £> increases. The MLE algorithm then will in essence 
"misuse" £, b c , and x Q to fit modeling error; this explains why generally 
better results were obtained when only the elements of ^ were identified. 
Thus, the acceleration contributions arising from inadequate trim conditions 
appear to be a more likely source of modeling error than model structure con- 
siderations associated with the rank of p or with the models for 'x^ and 
Yi assumed in equations (16) to (18). Evaluating the effect of such errors 
on the computed response time histories, as well as on the accuracy of <p, 
accordingly represents an important part of the overall model verification 
task. 


RESULTS AND DISCUSSION 

The principal steps in verifying the linear perturbation model of refer- 
ence 2 were to compare the respective estimated and measured output response 
time histories y^ and y mf i and to compare the corresponding analytical 
derivatives with the EKF and MLE values identified from the flight data. These 
comparisons accordingly involve three sets of y^ which were generated using 
the analytical derivatives of reference 2, the EKF values of reference 4, and 
the MLE estimates of $ from equation (2) . As stated in reference 4, the EKF 
derivatives were obtained using a six degree-of-f reedom model subsequent to pre- 
processing the flight data with a Kalman filter based on kinematical relation- 
ships. The effect of the SAS is taken into account in the EKF algorithm by 
adding the measured SAS signals directly to the control inputs in the equiva- 
lent open loop fashion described in connection with equations (10) and (11). 
Since the analysis being presented assumes decoupled three degree-of-f reedom 
dynamical modes, the three sets of y^ for each of the longitudinal and 
lateral-directional runs listed in table V were obtained separately by evaluat- 
ing equation (15) or (17) as appropriate. The corresponding time histories 
of x^ required in these computations were generated by numerically integrating 
equations (16) and (18) using the fourth-order Runge-Kutta method with a fixed 
step size of 0.10 sec. 

In order to avoid an unwieldy volume of plotted and tabulated data, which 
contributes little additional information, only results of runs from test 
points 2, 3, 5, 7, 10, 13, and 1 5 of table V are presented. An exception is 
run 33 from test point 11, which is included in connection with model structure 
verification. The results from the single input runs in test points 2, 3, and 5 
that involve ^coll an< 3 ^rud were also omitted as being unnecessary for the 
purposes herein. All of the manual input runs were eliminated from further con- 
sideration in the early stages of the derivative identification computations 
because these data did not yield convergent solutions for tp. Most of these 
runs resulted in y m ^ values that exceeded acceptable small perturbation 
bounds and generally exhibited poor excitation of both the short and long period 
response modes. 

In generating the MLE solutions, the reference 2 derivatives were always 
used for <p 0 in equation (2) since they represented the best available starting 
values for initiating the iteration process. Because the y m ^ responses to 
single input maneuvers contain no information about the control derivatives 
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corresponding to the omitted input other than through feedback (see eqs. (1 5) 
to (18)), these derivatives were not estimated and were held constant at their 
reference 2 initial values for all such runs. This fact is indicated by the 
parenthetically marked entries in tables II and IV. 


Longitudinal Derivatives and Response Time Histories 

The time histories y£ and y m £ for the seven longitudinal runs 5, 9, 

17, 25, 31 , 39, and 43 from test points 2, 3, 5, 7, 10, 13, and 15 of table V 
are plotted in figures 4(a) to 4(g), and the corresponding derivative values and 
trim conditions are listed in table II. No MLE results are presented for runs 9 
and 25 because convergence of tp was apparently prevented by the data problems 
discussed previously. Comparisons with the EKF method are made for only five of 
the seven runs since derivative values for the other two are not given in refer- 
ence 4. To investigate the validity of the model structure assumed in equa- 
tion (16) for the longitudinal mode, the results for run 33 of test point 11 
from two different flights are compared in figure 5 and table III. 

By referring to figure 4, reasonable agreement between the y£ generated 
from the analytical derivative values of reference 2 and the y m £ flight data 
is seen to exist for all seven runs, with some exceptions. These concern the 
presence of unmodeled response motion and discrepancies in the y£ response 
characteristics associated with the values of the modeled stability and control 
derivatives. Mention should be made of the fact that pilot model inputs as 
given by equation (22) were inadevertently added to 5-£ on g for runs 33 
(flight 035), 39, and 43. Other than to slightly distort the sinusoidal struc- 
ture of S lon g, the addition of ^longp££ ot to these runs did not appear to 

create problems or affect derivative identification. 

U nmod e led response .- The time histories plotted in figure 4 show an oscil- 
latory contribution to the response motion that occurs at nearly the damped 
natural frequency of the longitudinal short period mode and is not predicted 
by equation (16) for either of the three sets of derivatives listed in table II. 
This unmodeled response appears mainly in the n x plots but is also evident 
in those for q, n z , and Sq SA g for some runs. These effects, which are most 
prominent in runs 5, 17, and ^5, are almost absent from runs 31, 39, 43, and 
both runs 33 in figure 5. Although hysteresis and nonlinearity in the CH-47 
control system (or perhaps rotor dynamics) represent possible sources, the cause 
of the unmodeled response motion was not determined. 

Variation in longitudinal response characteristics .- The second area of 
disagreement exhibited by the y£ plots in figure 4 pertains to differences 
between the characteristics of the measured and computed responses to mainly 
the short period component A-j sin t of both 6£ on g and 5 CO £;£ inputs. 

(See eq. (21).) The results for all seven runs show the responses to 
Ai sin w-| t of each type input, as generated using the reference 2 derivatives, 
to occur at slightly higher frequencies and with smaller amplitudes than the 
measured data. Although there may be seme contribution to these differences 
from unmodeled accelerations arising from unsteady trim conditions, incorrect 
stability and control derivatives also can affect the response characteristics 
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of the computed time histories. The fact that the plots in figure 4 show dif- 
ferences between the from the reference 2 and MLE derivatives strongly sug- 

gests that error in these parameters probably accounts for most of the dis- 
agreement with the y m ,i data. 

In order to determine which derivatives dominate these variations in 
response characteristics, several families of y^ sensitivity curves were gen- 
erated in which each of the 15 derivatives modeled in equation (16) were set to 
zero one at a time. These plots, which were computed frcm the reference 2 deriv- 
atives for airspeeds of 0, 20, 40, and 60 knots, indicate that the frequency and 
damping characteristics of the y^ are most sensitive to changes in Mq. Refer- 
ence to table II suggests that the reference 2 values for Mq are too large, 
which thus appears to be the main reason the frequency and damping of the corre- 
sponding short period responses are accordingly too large. The generated 

from the MLE derivatives tend to agree more closely with the y m but the 
values for Mg may be somewhat too small, particularly at the higher airspeeds. 
Except for runs 17 and 25, the EKF results also indicate the reference 2 values 
for Mq are too large and agree better with the MLE estimates. 


In addition to showing the short period damping to be dominated mainly by 
Mq, the sensitivity plots further indicate the amplitudes of the y^ responses 


'■'long and ‘-'coll inputs to be approximately proportional to the magnitudes 


to 6 long and 

of the corresponding control derivatives, as expected. The pitch rate time his- 
tories plotted in figure 4 thus indicate that the reference 2 and MLE values 


for M<$ 


long 


listed in table II are of about the correct magnitudes, while those 


for the EKF model appear to be three to four times too small. Similarly, the 
plots of n x suggest that the MLE values for x 6 long are reasonably correct 


and that those for the other two models are an order of magnitude too small. 
The plots of n z exhibit about the same behavior for the runs at hover and 
20 knots, except that the sign of z S]_ on also appears to be incorrect for 


the reference 2 and EKF models. 


The three sets of 


Z ^long 


values are in better 


agreement for the runs at 40 and 60 knots; however, the sensitivity plots indi- 
cate that damping of the n z responses is also strongly affected by Z w . These 
results show the MLE estimates of Z w for runs 31 and 39 to be of incorrect 
sign and imply that the reference 2 values are perhaps a factor of 2 too large. 
The 6 qqh control derivatives for the three derivative models agree more 
closely than those for ‘Siong and do not give rise to any large differences 


between the amplitudes of y^ and y m , i* 


Pitch SAS model .- The plots for &Q r SAS i n figure 4 (and also in fig. 5) 
generally exhibit the same response characteristics as would be expected for q 
since the SAS states v-\ , V 2 , and V 3 are dominated by the pitch rate response. 
(See eqs. (12) and (16).) These curves indicate about the same level of agree- 
ment between the measured pitch SAS time histories and those predicted by each 
of the three sets of derivatives. Thus, the <50 f SAS results tend to verify 
the pitch SAS modeling employed in equations (15) and (16). 
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Model structure .- The comparison of results from two different trials of 
run 33 from test point 11, presented in figure 5 and table III, provides data 
for investigating the validity of the model structure assumed in equation (16) 
and represents an important aspect of the overall model verification task as 
previously discussed. In figure 5(a), and vice versa in figure 5(b), the y^ 
generated by the $ identified using the plotted y m ^ are compared with both 
the y^ computed with the $ from the other run and the y^ from the refer- 
ence 2 derivatives. Both sets of results indicate relatively minor degradation 
of the curve fit when $ from a different trial of the same run is used to gen- 
erate y^. The agreement between y m •[ and the MLE estimates of y^ from both 
runs indicated in figure 5 and table III does not, therefore, appear to be just 
the result of satisfying the MLE minimization criteria but represents limited 
verification of equation (16). As the trim conditions for both runs were gener- 
ally better than those for most of the other higher airspeed runs, the degrada- 
tion in accuracy of the $ due to nonzero accelerations in the y m ^ is prob- 
ably not very great. 


Lateral-Directional Derivatives and Response Time Histories 

The time-history comparisons for the seven lateral runs 6, 10, 18, 26, 32, 
40, and 44 from test points 2, 3, 5, 7, 10, 13, and 15 are plotted in fig- 
ures 6(a) to 6(g), and the corresponding derivatives and trim conditions are 
listed in table IV. As with runs 9 and 25, data problems apparently prevented 
convergence of $ so that no MLE results for run 1 8 are presented. The EKF 
results for run 18 are also omitted since they appeared to be similarly affected 
by the same data problems. Comparisons with the EKF method accordingly were 
made for just four of the seven runs because the necessary derivative values 
for runs 1 0 and 44 are not given in reference 4. Unfortunately, none of the 
MLE solutions for the lateral runs provided suitable duplicate data for model 
structure verification of equation (18), as was done with equation (16). 

The results presented in figure 6 generally indicate about the same level 
of agreement between y^ from the reference 2 derivatives and the y m j_ flight 
data as for the longitudinal runs and exhibit similar discrepancies as well. 
Although no unintentional pilot model implementation occurred with the lateral 
runs, the fact should be mentioned that the S]_ at inputs for runs 26 and 32 
contain an unprogramed signal of undetermined origin which is not modeled by 
equation (21 ) . This signal resembles the pilot model input given by equa- 
tion (22); however, the amplitude variations in u and w over the relevant 
portions of their time histories are too small in both runs to generate a 
6iongp^ lot comparable to the input signal in question. While the source of 

this signal was therefore not the pilot stabilization model, its presence 
appeared to have no noticeable effect on the results. 

Unmode led response .- The time history plots in figure 6 show an unmodeled 
oscillatory contribution to y^ similar to that noted in figure 4; however, 
the frequency is generally higher relative to that of the lateral short period 
mode and the amplitude is smaller. As with the longitudinal runs, the unmodeled 
response is most prominent in the linear acceleration plots. This behavior, 
which is not predicted by equation (1 8) for either of the three sets of deriva- 
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tives listed in table IV, is evident in the ny curves for all seven runs and 
in the 6 ^ gAS Plots for runs 6 and 10 also. The period of this oscillation 
appears to be about 1 sec or roughly half that observed for the ones in the 
longitudinal runs. As in the case of the longitudinal mode, the cause of the 
unmodeled response was not identified. 

Variation in lateral-directional response characteristics.- The charac- 
teristics of the measured and computed responses to the A] sin 0 )-] t component 
of <$iat and 6 ru( j inputs exhibit differences that are similar to those noted 
for the longitudinal runs. In the same manner as for equation (16), sensitiv- 
ity time histories were generated for the 1 5 derivatives modeled in equa- 
tion (1 8) . These plots were also computed for airspeeds of 0, 20, 40, and 
60 knots using the reference 2 derivatives. They indicate that variation in 
Ip has the most effect on the lateral-directional short period frequency and 
damping characteristics. 

Examination of the plots in figure 6 and the entries in table IV sug- 

gests that the reference 2 values for Lp are perhaps 20 percent too small for 
the hover and 20 knot runs, and roughly the same amount too large for the runs 
at 40 and 60 knots. By comparison, the MLE estimates of Lp yield yj_ that 
agree more closely with the y m ^ at the lower airspeeds, but the entries for 
the higher airspeeds appear slightly too small except for run 44. With the 
exception of run 40, all of the EKF values for Lp are of positive sign as 
noted in reference 4. The positive values for L p result in insufficient damp- 
ing of the lateral-directional short period mode and cause some of the y^ in 
figure 6 to exhibit divergence. (See fig. 6 (a) for example.) 

With regard to amplitude, the principal areas of disagreement between y^ 
and y m , i concern the curves for p and n y plotted in figure 6 . The deriv- 
ative sensitivity plots indicate that the differences in the response amplitudes 
of p to inputs are not due to the damping action of Lp alone. The 

plotted results thus imply that the values for L$ listed in table IV are 

X at 

too small by 30 percent or more for all three derivative models. The curves 
for n y show 180° phase differences in the responses to both 6 ^a t and ^rud 
inputs that indicate the reference 2 and EKF values for the lateral force con- 
trol derivatives have incorrect signs. The entries in table IV and the n y 
plots generated from the MLE derivatives suggest that Y$ should be nega- 

tive for all seven runs with approximately the magnitudes of the reference 2 
values. Similarly, y S rU( 3 should remain positive with magnitudes somewhat 

larger than the MLE estimates. 

Roll SAS and yaw SA5 models .- In a manner analogous to the 69 SAS plots 
in figures 4 and 5, the time histories of <5^ SAS and 6 ^ SAS presented in 
figure 6 exhibit about the same response characteristics as those for p and 
for the composite of r and n y , as would be inferred from equations (13), 

(1 4) , and (1 8) . The measured and computed SAS signals again reflect the same 
level of agreement as that between the y m ^ and the three sets of y^, except 
that the amplitudes of the roll SAS measurements are 10 to 20 percent larger 
than predicted by either of the three derivative models. The reason appears 
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to be that canceling one of the two roll SAS actuators to halve the roll SAS 
gain Kp (see table I) , as recommended in reference 3 to minimize rotor- 
fuselage coupling predicted by simulation studies, did not yield exactly the 
desired 6.35 value. However, the results plotted in figure 6 do not otherwise 
indicate any deficiencies which would invalidate the SAS models incorporated 
in equations (17) and (18). 


Closed Loop MLE Derivative Identification 

As mentioned in the Introduction, one of the objectives of the overall 
model verification was to investigate the effect on derivative identification 
of conducting the flight tests on a closed loop basis with the SAS operating. 

The preliminary checkout of the MLE algorithm for the CH-47 application, using 
simulated y m ,i data which include the effects of the SAS, indicated that good 
closed loop identification accuracy should be attainable with reasonably consis- 
tent flight data. The M LE results presented in figures 4 to 6 and tables II 
to IV strongly suggest that comparable accuracy also would have been achieved 
with the actual flight data if the trim conditions had been more adequately 
obtained. Although the modeling error introduced by inadequate trim conditions 
may retard or even prevent convergence of the MLE algorithm, the checkout compu- 
tations further showed that no increase in the number of iterations of equa- 
tion (2) should occur because of performing derivative identification with a 
closed loop dynamical system. 

The notions that the SAS will suppress response modes and cause what is 
termed masking of the derivatives, so as to reduce or impair their identifi- 
ability, were not observed to be true of the $ computations with either simu- 
lated or actual flight data. To the contrary, the plots of y m ^ in figures 4 
to 6 show well defined responses with no evidence of any suppression of either 
the short or long period modes. By preventing divergent responses, the SAS thus 
appears to facilitate rather than impede good modal excitation. In addition 
to decoupling the longitudinal and lateral-directional modes very effectively, 
as discussed in conjunction with equations (15) to (18) , the SAS also seems to 
enhance the identif iability of the control derivatives when the MLE algorithm 
is used. 


Comparison of Manual and Computer-Generated Input Maneuvers 

The primary purpose in comparing results from manual and computer-generated 
input maneuvers is to evaluate the suitability of each for generating flight 
data for derivative identification purposes. In this connection, the most 
important requirement for the response data is good excitation of the aircraft 
dynamic modes at the largest amplitudes that remain within practical small per- 
turbation bounds. The y mf i curves for run 36 plotted in figure 7 and for 
run 40 in figure 6 illustrate typical response measurements generated by manual 
and computer-generated inputs at the same test conditions. Although the plots 
in figure 7 are not representative of typical pilot performance in any sense 
because no practice maneuvers were flown, they do however demonstrate inherent 
difficulties that can be encountered in manually implementing a doublet pulse 
input maneuver having the desired frequency and amplitude characteristics. In 
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this example, the amplitude of the 6i at input is satisfactory but the fre- 
quency is almost three times less than the value for listed in table VI. 

Primarily because the pulse duration is thereby too long, the resulting ampli- 
tude of p (see fig. 7) is too large and is seen to be roughly double that 
produced by the corresponding computer-generated 5i a t input maneuver. The 
manual inputs also are not very repeatable, as both the pulse duration and 
amplitude are usually determined by the pilot's judgment with the result that 
accurate repetition of even the most simple doublet pulse is almost impossible 
to achieve. Furthermore, the complexity of the piloting task essentially pro- 
hibits using inputs that require superimposing multiple pulses of differing 
frequency and amplitude content. 

Computer-generated inputs, by comparison, exhibit the potential for alle- 
viating practically all of the difficulties just described. The results pre- 
sented in figures 4 to 6 indicate that these inputs can solve the problem of 
response mode excitation very successfully. Most of the y m ^ curves show 
well-defined sinusoids that have precisely the intended short and long period 
response frequencies. In this connection, using the variable scale factor K 
(see eq. (22)) proved to be highly effective in constraining the response ampli- 
tudes of the Ym,i in these runs within practical linearity limits. 

The onboard computer, by providing the capability for programing any vari- 
ety and number of input shapes or waveforms that can be represented mathemati- 
cally, thereby offers the means for assuring proper modal excitation which effec- 
tively solves the most difficult part of the input design problem. The remaining 
task then concerns determining an optimum combination of the number of sinusoids 
or other input shapes, their frequencies and amplitudes, and the number of 
cycles each one is to be repeated during the course of the input maneuver. Since 
they can be precisely repeated as often as desired once programed, these inputs 
also can facilitate generating duplicate y m ,i data for model structure verifi- 
cation. Thus, when the test aircraft is equipped with the necessary electronic 
control system, the use of computer-generated inputs appears to be superior to 
manual ones. Subsequent flight tests for other purposes showed that very good 
trim conditions can be achieved by means of the electronic control system. 

Hence, the problem with nonzero trim accelerations can be largely avoided by 
also automating the trimming of the aircraft prior to initiating a test 
maneuver . 


CONCLUSIONS 

The results from flight tests for verifying an existing analytical model 
of the aerodynamic derivatives for a CH-47 tandem-rotor transport helicopter 
indicate the following conclusions: 

1 . The analytical model is reasonably correct with respect to the stability 
derivatives, which exhibit only minor inaccuracies, but some of the longitudinal 
and lateral force control derivatives have incorrect signs and magnitude errors 
as large as a factor of 10. 

2. Although nonzero trim accelerations in the flight data caused the deriv- 
atives identified from both extended Kalman filter (EKF) and maximum likelihood 
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estimator (MLE) algorithms to be in error also, the results from the latter 
show a strong tendency to correct the deficiencies in the analytical model that 
indicates the MLE algorithm would yield good accuracy if trim accelerations were 
absent. 


3. Closed loop flight testing with the aircraft stability augmentation sys- 
tem (SAS) operating does not adversely affect derivative identification for 
tandem-rotor helicopters. 

4. The use of computer-generated input maneuvers is superior to manual 
ones and is highly effective in generating transient response data having good 
modal excitation without exceeding practical small perturbation amplitude 
bounds . 


5. Whenever the necessary electronic control system is available, the prob- 
lem with nonzero trim accelerations should be largely avoidable by also automat- 
ing the trimming of the aircraft prior to initiating a test maneuver. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
November 21 , 1 979 
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TABLE I.- CH-47 SAS MODEL CONSTANTS 


Feedback gains: 

Kg, cm/r ad/sec 

a Kp, cra/r ad/sec 

K r , cm/rad/sec 

K r,p' cm/rad/sec 

K V £ 

U 0 = 0 to 20 knots, cm/m/sec 
U 0 = 20 to 40 knots, cm/m/sec 

U D = 40 to 60 knots , cm/m/sec 

U Q = 60 to 80 knots, cm/m/sec 


State variable coefficients: 


b 


a 0 

a 1 

a 2 

a 4 


a 6 


Feedback gain coefficients: 

b 1 

b 2 

b 3 

b 5 

d 7 


43.94 

6.35 

26.21 

16.43 

0 

0.1233 

0.2717 

0.2958 


0.0364 

0.8177 

3.0401 

0.3125 

0.2433 


-1 .2890 
3.2676 
-8.8799 
0.3125 
0.2433 


a Value one-half that given in reference 5. (See also ref. 3.) 
b a4 = 2a4; U Q < 40 knots. 
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TABLE II.- LONGITUDINAL AERODYNAMIC DERIVATIVES AND TRIM CONDITIONS 


Quantity 


U Q , knots 
W Q , ra s _1 


U 0 , knots 

W. . m 0-1 


deg 


Stability: 

Mq, S' . 
Mg, s' . 

Mq f s" 1 . 


M w , 

M*, S' 
M w , S' 

Mu. 

M u , 

M„, 


r 1 m - 1 

-1 m- 1 
1 nr! 

s _1 m -1 
s" 1 

. 0-1 m -1 


z w , 

^Wf 


Xy' 

x w , 

x w , 

x u , 

Xu r 
X,„ 


s 1 
s-1 
s ^ 

s -1 

s -1 

S~1 

S _1 

s“l 

o-l 


Control: 
Mfi 


ms: 

MS 


long' 
long ' ^_2 
long' 


cm 

2 _-l 


Mfi 

MS, 

MS 


coll' 
ooii' s _: 


ooll ' 


cm 

2 


s -2 an 1 
=-2 


^ cm 1 
2 


Derivative 

model 

(a) 


O 

□ 

o 

o 

□ 

o 

o 

□ 

o 

o 

□ 

o 

o 

□ 

o 

o 

□ 

o 


o 

□ 

o 

o 

□ 

o 

o 

□ 

o 

o 

□ 

o 

o 

□ 

o 

o 

□ 

o 


1.7 

-o.oo 

6.36 


-1 . 2293 
-0.8882 
-0.6680 

0.0077 

-0.0284 

-0.0496 

0.0304 

0.0321 

0.0344 

-0.2956 

0.1203 

0.0454 

0.0326 

0.2425 

0.0382 

- 0.0211 

-0.0690 

-0.0429 


0.1296 
0.1116 
0.0588 

0.0075 
k (0. 0075) 
-0.0055 

0.0036 
-0.1608 
0. 0037 

-0.9674 
b (-0 . 9674) 
-0.7836 

0.01 37 
0. 1 51 4 
0.01 41 

0.1126 
b (0. 1126) 
0. 0980 


0 

2.54 


1.0 
2.1 5 
7.42 


-1 .2486 


0.0094 


0.0289 


-0.9661 


0.0131 


0.1125 


a O Reference 2 □ MI£ 

“Held constant at reference 2 values. 


O EKP (ref. i) . 


Run 1 7 


20 

2.54 


1 8.4 
3.67 
3.90 


-1.3422 

-0.6005 

-0.1471 

0.0402 
0. 2248 
-0. 0483 

0. 031 7 
0. 01 26 
0. 0009 

-0. 31 55 
-0 . 41 21 
-0.5977 

0.0322 
0.1 947 
0 . 01 80 

0.0013 
-0.01 26 
-0.0323 


0.1289 
0.1577 
0.01 29 

0.001 3 
b (0.0013) 
-0.0033 

0 . 01 99 
-0.1466 
0. 01 42 

-0.9784 
b (-0 . 9784) 
-1.1295 

0. 01 26 
0.1 784 
0. 0269 

0. 0831 
b (0. 0831) 
0. 0743 


Run 25 

40 

0 


44.9 
- 0.00 
3. 20 

-1.5695 
-1 .6089 
0.0583 
-0. 0363 
-0.0034 
0. 001 4 
-0.4454 
-0.0429 
0.0387 
0.031 5 
- 0.011 6 
-0.0358 

0.1 401 
0.0206 
0.0359 
0.011 5 
0. 0444 
0.0206 
-7.0282 
-1 .1296 
0.01 30 
0.01 43 
0. 0655 
0.0655 


Run 31 
40 

2.54 


41 

.2 

2 

!.78 

5 

i. 31 

-1 

.6084 

-0 

.4776 

-0 

.4633 

0 

.0680 

-0 

.0591 

-0 

.01 57 

-0 

.0076 

-0 

. 0239 

-0 

. 0008 

-0. 

.41 55 

0, 

.41 07 

-0. 

.1 055 

0. 

.0385 

-0. 

,011 7 

0. 

0436 

-0. 

0093 

-0. 

0383 

-0. 

01 54 

0. 

1 395 

0. 

1 080 

0. 

031 3 

0. 

0392 

0. 

0137 

0. 

0254 

0. 

0512 

0. 

1184 

0. 

0604 

-1 . 

001 3 


-0.7883 

-0.6787 

0.01 25 
0.1363 
0.01 21 

0.0661 

0.0444 

0.0642 


Run 39 


60 

0 


62.8 

0.93 

2.49 


-1 .681 8 
-0.51 59 
-0.6040 

0.0579 
0.0480 
-0.01 07 

-0 . 01 38 
-0.001 3 
-0 . 01 57 

-0.551 2 
0.0758 
-0.1 550 

0.0376 

0.0755 

0.0201 

-0.0205 

-0.0304 

-0.0483 


0.1 540 
0.1 397 
0.0473 

0.0602 
0.01 87 
0.011 7 

0.0561 
0.11 06 
0.0853 

-1 .1232 
-0.6667 
-0.781 2 

0.01 52 
0.1 744 
0.01 37 

0.051 2 
0.0568 
0.0523 


Run 43 


60 

2.54 


63.2 
3.71 
1 .91 


-1.7072 

-0.2911 


0.0623 
0. 0097 


-0 . 01 70 
0. 01 00 


-0.5345 

-0.2266 


0.0390 
-0. 0340 


-0. 01 92 
0.0001 


0.1 549 
0.1 278 


0. 0662 
0.0306 


0.0630 

0.0355 


-1.0990 

-0.7668 


0. 01 50 
0.1257 


0.0550 

0.0670 


TABLE III.- LONGITUDINAL AERODYNAMIC DERIVATIVES AND TRIM CONDITIONS 
FOR RUN 33 FROM TWO DIFFERENT FLIGHTS 


Quantity 


knots 
« 0-1 


Jcoli’ 

'-’long' 

^coll' 

^long' 

X(S coll' 


Reference 2 

o 


40 

5.08 

4.68 


-1 .6547 
0.0790 
-0.01 49 
-0.3849 
0. 0384 
-0.0069 


0.1 397 
0. 0450 
0.0602 
-0.9680 
0.011 9 
0. 0663 


Flight 031 

Flight 035 

□ 

o 

41 .7 

39.6 

5.52 

6.61 

3.55 

3.66 

-1 .1 070 

-0.5236 

-0.11 95 

-0.1 030 

0.1 803 

-0.0262 

-0.361 2 

-0.4440 

-0.1 071 

-0 . 1 21 7 

0.21 45 

-0. 0305 

0.0911 

0.0890 

0. 01 02 

-0.0035 

-0.1 61 5 

-0.1 548 

-0.8264 

-0.7306 

0.1 308 

0.11 55 

0. 0354 

0.0555 


29 



TABLE IV.- LATERAL-DIRECTIONAL AERODYNAMIC DERIVATIVES AND TRIM CONDITIONS 


Quantity 


















Run 6 

Run 1 0 

Run 1 8 

Run 26 

Run 32 

Run 40 

Run 44 






(a) 








0 o , knots 






0 

0 

20 

40 

40 

60 

60 

W Q , m s* 1 






0 

2.54 

2.54 

0 

2.54 

0 

2.54 

Trim: 













U Q , knots . . . . 






0.4 

1 .7 

21 .8 

38.4 

47.1 

62.8 

63. 4 

W Q , m s _1 . . . 






-0.00 

2.48 

3.20 

-1 .71 

2.33 

0.85 

3.51 

®o» deg 






8.48 

6.72 

4.72 

4.25 

3.08 

3.20 

0.62 

Stability: 













*V. S-] 





O 

-0.1 371 

-0.1 31 8 

-0.1 083 

-0.0689 

-0.0642 

-0.0740 

-0.0704 







-0.11 70 

0.1 531 






1 . 

O 

o 


0.1872 
-0.01 82 

0.1393 
-0.01 71 

-0 . 1 01 0 
-0.01 80 

U.UU4I 

0.0481 

-0.0200 


L v , s 1 m 1 . . . 





-0. 01 91 

-0.01 85 

-0.01 68 






□ 

o 

-0.0202 

0.0216 

-0.0851 








-0.0239 

0.0056 

0.0051 


v ' 








V s "! 





o 

-0.6950 

-0.731 3 

-0.7489 

-0.7748 

-0.8270 

-0.81 84 

-0.8789 







-1.3592 

-1 .1 532 




-0.6662 

-0.9739 















O 

0. 4069 



0.9993 

0. 2083 

-0.7586 










N v , s _1 m 1 . . . 





o 

□ 

o 

-0.0020 
-0.0000 
-0.01 58 

-0.001 8 
0. 0076 

0.0043 

0.0009 

0.0022 

0.0023 

0.0007 
-0 . 011 7 
0.0007 

-0.0003 

0.0082 

-0.0080 

-0.0005 

-0.0040 

... 










Hp- s _1 





o 

0.0008 

-0.001 5 

-0.0098 

-0.01 44 

-0.01 42 

-0.01 66 

-0.01 50 






□ 

-0.2458 







"-i 



-0.0898 

-0.0344 

0.0396 

‘ 

O 

-0.0301 

P' “ 








s-1 





o 

-0.0417 

-0. 0408 

-0.0409 

-0.0398 

-0.0394 

-0.0391 

-0. 0383 







0. 4230 







v ■ ; 

o 



-0.0575 

-0.0472 

-0.0318 

. UO-i / 

-0. 0470 









Control: 













Y<$, . . m s“ 2 cm~^ 
Y«^ f m s-2 cm-1 
Y« . m s -2 cm T 





o 

0.1 391 

0.1 389 

0.1 371 

0.1 352 

0.1 343 

0.1 344 

0.1 331 





□ 

o 

-0.1 736 
-0.081 8 

0.0393 


-0.0531 
0.1 824 

-0.0477 
0.1 797 

-0.0924 
0.1 210 

-0.071 6 








Y W m s i Cm 'i 1 

Y «rua' m s : 2 
Y «rud, m 3 “ 





o 

-0.0065 

-0.0066 

-0.0045 

-0.0046 

-0.0050 

-0.0064 

-0.0072 





□ 

O 

b (-0.0065) 

-0.0059 

b (-0.0066) 


0.1 089 
-0.0036 

0.0792 

-0.0049 

0.0642 

-0.0063 

0.0509 








r L «lat S J =”_ • 

L «lat, 3 cm • 





o 

0.1 636 

0.1 634 

0.1 61 7 

0.1 601 

0.1594 

0.1 597 

0.1587 





□ 

O 

0. 2030 
0.1 637 

0.1 303 


0.1777 

0.2402 

0.1772 
0.1 635 

0. 2367 
0.0522 

0. 2065 








r L «r U a, s : 2 “i] ■ 

Lfi , s z cm 1 . 

^rud, _2 -T 

L6 ' s * cm . 

u rud. 





o 

-0.0547 

-0.0547 

-0.0523 

-0 . 051 8 

-0.051 8 

-0.0530 

-0.0532 





□ 

o 

b (-0.0547) 
-0.01 55 

b <-0.0547) 


-0.0558 

-0.0667 

-0.041 6 
-0.1 403 

-0.0489 

-0.0501 

-0.1 020 








N<5- s -2 cm' 1 . 

N6, 1 !"' s-2 cm-’ . 
N «la t ; s ' 2 Cm ' 1 • 





o 

0.0036 

0.0036 

0.0040 

0.0039 

0.0038 

0.0035 

0.0033 





□ 

o 

0.0204 

-0.0068 

0.031 8 


0.01 21 
-0.0036 

0. 01 29 
0.0052 

-0.0360 

-0.0007 

0.01 28 








N <5 s -2 cm -1 . 

°r ud , 2 <^-1 





o 

0. 0802 

0. 0801 

0.0792 

0.0780 

0.0775 

0.0776 

0. 0768 









0.0609 




„°rud, s _ 2 CIn _ 1 
Nfi , s L cm 1 . 

u rud. 





o 

0. 021 1 



0. 0490 

0.0291 

U.UGU 

0. 031 2 









a O Reference 2 □ MLE <> EKF (ref. 4). 

hHeld constant at reference 2 values. 
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TABLE V.- CH-47 FLIGHT TEST RUN SCHEDULE 


(a) Single input maneuvers 


Control 

input 


Airspeed U Q , 
knots 


Descent rate W D , 
m/sec 


Wind velocity, 
knots 


Altitude, 

m 


6 

6 

6 

6 


long 

lat 

rud 

coll 


^long 

5 lat 

^rud 

^coll 


6 

6 

6 

5 


long 

lat 

rud 

coll 


5 

6 
6 
6 


long 

lat 

rud 

coll 


6 

6 

6 

6 


long 

lat 

rud 

coll 


6 

5 

6 
6 


lat 

rud 

coi: 


Computer generated 


0 

(Hover) 

0 

3 

61 .0 

0 

2. 54 

3 

137.2 

0 

(Hover ) 

0 

4 

15.2 

(In ground effect) 

20 

2.54 

7 

259.1 


Manual 



0 

(Hover) 

0 

3 

61 . 0 

40 

0 

) 3 

243.8 


^lon< 

‘-’lat 

5 rud 

5 coi: 



TABLE V.- Concluded 


(b) Computer-generated double input maneuvers 


Control 

Airspeed U Q , 

Descent rate W Q , 

Wind velocity, 

input 

knots 

m/sec 

knots 

ng? ^coll 

40 

0 

1 3 

t ? ‘-’rud 



1 0 

ng? ‘-’coll 

40 

CO 

O 

in 

1 

8 

t? 5 rud 



7 

ng? ^coll 

40 

-2. 54 

8 

t ? 6 rud 



8 

ng? ‘-’coll 

40 

2. 54 

11 

t? 5 rud 



11 

ng? ^coll 

40 

5.08 

11 

t ? 5 rud 



8 

ng? ^coll 

60 

0 

20 


20 

6 

6 


Altitude 

m 

243.8 
21 3.4 
228.6 
259.1 
274.3 
243.8 


20 

8 








On On On 


TABLE VI.- COEFFICIENTS FOR COMPUTER-GENERATED INPUTS 



Time, sec 


u 


a6 


long 
coll 
lat 
r ud 


Al , 

A 2 r 

0>J , 

w 2 , 

to, 1 ' 

to, 2 r 

ti , 

t 2r 

cm 

cm 

rad/sec 

rad/sec 

sec 

sec 

sec 

sec 


Single input maneuvers at hover and 20 knots 


— 

5. 08 

1 .27 

3.43 

— 

0.53 

0 

0 

1 .83 

1 2 

3.81 

1 .27 

3.43 

.53 

0 

0 

1 .83 

1 2 

3.81 

. 25 

2.70 

.32 

0 

0 

2. 33 

>1 5 

3.56 

.51 

2.70 

.32 

0 

0 

2.33 

>1 5 


Double input maneuvers at 40 and 60 knots 


6 

6 

6 

5 


long 
coll 
lat 
r ud 


3.56 

0.89 

3.43 

0. 53 

0 

6 

1 .8 

1 2 

3.56 

.89 

3. 43 

.53 

6 

6 

7.8 

1 2 

2.54 

.25 

3.50 

.32 

0 

0 

1 .8 

>1 5 

-2.54 

-.25 

3.50 

.32 

6 

6 

7.8 

>1 5 


a Longitudinal inputs were implemented with opposite signs in flight 


test because of sign convention employed in reference 3. 
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Figure 1.- Digital filtering of flight data from run 10 




u> 

<y» 



Time, sec Time, sec 


Figure 3.- Typical trim conditions. Run 18; U c = 20 knots; W Q = 2.54 m/sec. 



2 


+ Measured 
o Reference 2 



1 1 1 1 1 n 1 1 1 1 1 1 u 1 1 1 1 il i it 1 1 ii 1 1 1 1 1 it 1 1 mil i ii ii ml 1 1 1 1 j 1 1 1 1 1 1 1 1 1 1 mi I m 1 1 1 l 1 1 1 1 1 ■ ■ ■ . ■ . m I 

0123456789 10 

Time, sec 



Time, sec 


(a) Run 5; U 0 = 0; W 0 = 0 . 

Figure 4.- Comparison of longitudinal response time histories. 





, g units 



Time, sec 


(d) Run 25; U 0 = 40 knots; W Q = 0. 
Figure 4.- Continued. 


40 



g units q, deg/sec 





012345 6789 10 01234 567 89 10 

Time, sec Time, sec 

(a) Run 33 (flight 031); U Q = 40 knots; W Q = 5.08 m/sec. 

Figure 5.- Comparison of longitudinal response time histories for 
two different flights at same test conditions. 









, deg/sec p, deg/sec 






g units 



-.12 [- 
-. 16 * 

cT“T 2 3 4 5 6 7 8 9 10 

Time, sec 

(c) Run 18; U Q = 20 knots; W 0 = 2.54 m/sec. 
Figure 6.- Continued. 
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+ Measured 
O Reference 2 
□ MLE 
5 EKFfref. 4) 




1 1 1 1 1 1 1 1 1 1 1 i I 


I I ij I 11 III llullI'UIIII I I I 1 1 L I I j I ll I I 1 I I lii'll I 1 1 1 f ■ I I { I I II I J 1 1 I I III 

456789 10 012 

Time, sec 

(e) Run 32; U 0 = 40 knots; W 0 


1 1 1 l 1 1 1 n 1 1 1 nlnm mill I 

3 4 5 6 

Time, sec 

= 2.54 m/s ec 


Figure 6.- Continued 





i ii 1 1 1] n unliiiiiiiiiliiumnlimiiinliiinuiilimniuJ iiiiiimilniiii 1 1 1 I iiJiiimlimiimliiiimiilj iiimiili 


4 5 6 

Time, sec 

7 

8 

9 10 0 

1 2 

3 4 5 6 

Time, sec 

(f) 

Run 

40; 

U 0 = 60 

knots; 

o 

II 

lE^ 


Figure 6.- Continued 





S units 0, deg r, deg/sec p, deg/sec 



Figure 6.- Concluded 
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